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Abstract 

We study the low transverse momentum (pr) distribution of the Z-boson at hadron colliders for 
Pt ~ Aqcd using a factorization and resummation formula derived in the Soft Collinear Effective 
Theory (SCET). In the region px ~ Aqcd-, new non-perturbative effects arise that cannot be 
entirely captured by the standard parton distribution functions, and require an additional new 
non-perturbative transverse momentum function (TMF). The TMF is field-theoretically defined 
in SCET, fully gauge invariant, and captures the non-perturbative dynamics that affects the px- 
distribution in the region px ~ Aqcd- The TMF also reduces to the expected perturbative result 
in the region px ^> Aqcd- We develop phenomenological models for these TMFs in the non- 
perturbative region and present example fits to the available data. 
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I. INTRODUCTION 



The description of the low transverse momentum (p?) distribution of electroweak gauge 
bosons and the Higgs boson has been the subject of extensive study |I|-|16|- It plays an 



important role in the precision measurement of the P^-boson mass and Higgs boson searches 
while providing an important test of perturbative Quantum Chromodynamics (QCD). In the 
region of low transverse momentum px *C M, where M denotes the mass of the electroweak 
gauge or the Higgs boson, large logarithms of pt/M spoil the perturbative expansion in the 
strong coupling and require resummation. 

More recently, the low-p^ resummation was studied using a factorization theorem de- 
rived [T7], |TB[ in the Soft-Collinear Effective Theory (SCET) [|T^-^T|. The result derived using 
SCET can be written entirely in momentum space, avoiding issues arising with the impact- 
parameter space present in the standard approach. All objects in the factorization theorem 
have well-defined operator expressions. A detailed study of the region Aqcd <C pr <C M 
was performed for the production of Higgs and electroweak gauge bosons. In this region, the 
factorization theorem is given entirely in terms of perturbatively calculable functions and 
the standard initial state parton distribution functions (PDFs), and takes the schematic 
form 

" ' a H®g®f®f. (l) 



dp T dY 



Convolutions between the various objects are denoted by the symbol <S>, H denotes a hard 
function whose renormalization group (RG) evolution sums logarithms of pt/M, Q denotes a 
perturbative function at the pr-scale and describes the emission of soft and collinear partons 
that recoil against the heavy boson, and / <8> / denotes the product of the initial state PDFs 
which are evaluated at the p^-scale as determined by DGLAP evolution. Resummation of 



the large logarithms was performed at the next-to-leading log (NLL) accuracy in Ref. [18 



using renormalization group (RG) evolution in the effective theory. The results for the Z- 



boson are in excellent agreement with Tevatron data collected by the CDF |22| and DO f23 
collaborations. 

In this paper, we turn our focus to the region p? ~ Aqcd <^ M where new non- 
perturbative effects arise that cannot be captured entirely by the standard PDFs. The 
region pt ~ Aqcd is sensitive to the transverse momentum distributions of the partons in 
the initial state hadrons and to transverse momentum emissions of order Aqcd- in this 
region, the factorization formula takes the schematic form 

d 2 a 

H®1C, (2) 



dp\ dY 



where the function K is evaluated at the scale \it ~ Pt ~ Aq C d- The definition of /C is given 
in Section [TTT| . It is a new non-perturbative function that cannot be described in terms of the 
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standard PDFs alone. In order to facilitate a smooth transition between Eqs. (0) and (P as 
one increases pt from non-perturbative to larger perturbative values, it is useful to write !C 
in the form 

ic~g®f®f, (3) 

which defines the new transverse momentum function (TMF) Q. For phenomenological 
purposes, the TMF is modeled in the non-perturbative-pr region with the constraint that it 
reduce to the expected perturbative result in Eq. (p]) in the high-p-r region. For this reason 
we use the same symbol Q to denote the TMF over the entire pt spectrum. The function 
tC, or equivalently the TMF Q, is universal and depends only on the hadronic initial state. 
The region of pt ~ ^qcd has been studied extensively in the context of semi- inclusive 



deep- inelastic scattering (SIDIS) [p4| - |28fl , and also within the Collins-Soper-Sterman (CSS) 
approach to resummation of low-pr logarithms [p9H3~lj. In SIDIS processes, transverse 
momentum dependent parton distribution functions (TMDPDFs) typically arise in order 
to describe the order Aqcd dynamics in the initial hadrons. The TMDPDFs are typically 
not gauge invariant under singular gauge transformations, and arriving at a gauge-invariant 



definition has been the subject of much research |24-f28|, |32|-j35|l. In our formalism, it is 
instead JC that is the fundamental non-perturbative object in the region pj< ~ Aqc£>. It is 
fully gauge invariant |I7j. We choose to write fC in the form of Eq. ([$]), and view the TMF 
(Q) and the PDFs (/) as the fundamental objects of interest. Both of these are manifestly 
gauge invariant, and have a more intuitive and smoother connection with the form of the 
factorization theorem in the region pt ^> Aqcd- Thus, instead of TMDPDFs it is simpler 
to work with the gauge invariant TMFs and PDFs which are well-defined over the entire 
PT-spectrum including both perturbative and non-perturbative values of pt- 

The goal of this manuscript is to develop initial models of the TMF that satisfy the 
following criteria: they reduce smoothly to the perturbative result as one increases pt, and 
they preserve the RG running of JC in order to cancel the running of the hard function, as 
required by the scale invariance of the cross section. We do so, and present numerical results 
for pp initial states by fitting to Tevatron data. A global fit to all available data and an 
analysis of fixed-target data in order to describe the pp initial state and test the universality 
of the TMF is reserved for future work. 

The outline of our paper is as follows. In Section O, we review the factorization formu- 



las for the perturbative px region derived in [O, [18] and briefly discuss the various pieces 



and notation. In Section |TTT| we present the factorization formula for the non-perturbative 
p T -region and discuss the issues involved in developing a non-perturbative model. We give 
phenomenological models for the TMF Q in the non-perturbative px region and show nu- 
merical results in Section [TV]. We conclude in Section |V[ 
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II. THE PERTURB ATIVE p T REGION 



In this section we briefly review the basic elements of the factorization and resummation 
formula for the transverse momentum distribution of the Z-boson in the region Aqcd "C 
Pt <C M, as derived in Refs. |17|, |18| . Although we focus on the Z-boson, the analysis is 
similar for any color-neutral heavy final state. The appropriate effective field theory for this 
observable is SCETn, which has both collinear and soft degrees of freedom that can recoil 
against the Z-boson with transverse momenta of order p?. The collinear and soft degrees 
of freedom have momentum scalings 

p n ~ M(r] 2 , l,r/), p n ~ M(l,r] 2 ,r]), p s ~ M(r),r),T)), rj ~ (4) 

iviz 

where we have used the notation p = (n • p,n ■ p,p±) to denote the light-cone and 
transverse momentum components. The light-cone four-vectors are = (1, 0, 0, 1) and 
n M = (1,0,0,-1). The p n ^ momenta denote collinear momenta with large components 
along the n M and n M directions respectively. The soft momenta are denoted by p s . The 
transverse momentum distribution in the region Aqcd <C Pt *C M is dominated by these 
collinear and soft modes radiated from the initial state partons. In SCETn, these emissions 
build up into collinear and soft Wilson lines that dress the Z-production current. The final 
factorization and resummation formula for the differential cross-section of the Z-boson as a 
function of its transverse momentum and rapidity (Y) is given by 



d 2 a 7T 2 r 1 r 1 dx r l dx[ r l dx' 2 

dp 2 T dY A c 2 J Xl J %2 J Xi x[ J X2 x' 2 



x HKx^Q 2 , hq\ /j t ) g qrs (x u x 2 , x[, x' 2 ,p T , Y, fi T )f r (x[,i2 T )f s (x' 2 , /i T ). 

(5) 

The above formula involves a convolution of three types of objects: the hard function H q z , 
the TMF Q qrs , and the initial state PDFs / r)S . The indices r, s run over the initial partons 
and the superscript q denotes the fact that the Z-boson production vertex involves a quark 
current. The hard function describes the physics of modes with virtuality p 2 ~ M 2 that 
are integrated out at the scale \Lq ~ M. The hard function is then evolved down to the 
scale fiT ~ Pt via its renormalization group equations, summing large logarithms of order 
M/pt in the process. The TMF function Q qrs lives at the ht ~ Pt scale and describes the 
physics of the soft and collinear emissions in a way that is consistent with the constraints 
imposed on the px and Y of the Z-boson. The initial state PDFs f r>s are evaluated at the 
/it scale after DGLAP evolution from the non-perturbative scale, summing logarithms of 
order Aqcd/pt in the process. 
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The TMF function Q qrs has the form 

d 2 b 



/d b f 
j^^Jo[b ± p T ] J dt+dt^l n . tqr (-±,t+,b ± ,fi T ) 

W^T. ««, b±, ^t)S- 1 (x 1 Q - e Y Jp 2 T + M 2 - \, x 2 Q - e~ Y J^ + M 2 - 



(6) 



where the functions T n . qr ,T n .q s correspond to collinear emissions in the n and n directions 
respectively and S correspond to soft emissions. The inverse soft function (iSF) <S _1 arises 



due to zero-bin subtractions |L7|, [18|, |36|-p^] necessary to avoid the problem of double-counting 
the soft region. The collinear function I n - qr is defined through the matching of a nucleon 
matrix element called the impact-parameter beam function (iBF) B% onto the standard 
PDFs as 

~ f dz x 

B*(x,t,b±,HT)= —Zn;qr(-,t,b ± ,H T )f r (z,H T ), (7) 

J x Z Z 

with an analogous equation for the n-sector. For precise field-theoretic definitions of the 
iBFs, the iSF, and the hard-function H q z we refer the reader to Refs. [ 17| , PI . Analogous 



nucleon beam functions p0|-[43| are known to appear in other collider processes. We note that 
this SCET formalism accomplishes the resummation of large logarithms differently than the 
traditional QCD approach. Logarithms of the matching- scale ratio ln(fj,Q / n?) that appear in 
the partonic cross section are resummed via the RG evolution of Hq. Upon identifying fiQ ~ 
M and fix ~ Pt, these become the standard small-p^ logarithms. Kinematic logarithms 
which directly have ln(M/p^) appear after integration over the momentum fractions x 1 2 in 
Eq. (^|). It was shown in Ref. that this formalism reproduces the correct logarithms upon 
expansion of the resummed result to the fixed order 0(a 2 ) given our current knowledge of 
Q in perturbation theory. 



We comment briefly on the recent work of Ref. [H] which also uses SCET to address 
Pr-resummation. We disagree with several aspects of their results. Their analysis is based 
on the claim that the emission of soft radiation with transverse momentum of order pt does 
not affect the spectrum of the Z-boson. This is in contrast to our effective field theory 
(EFT) where both collinear and soft radiation, with transverse momentum of order p?, play 
a dynamical role in determining the transverse momentum spectrum. It is well-known [|Hj|- 
U that the emission of multiple collinear and soft partons from the initial-state collinear 



partons build into eikonal Wilson lines and is a leading order effect in SCETn- Since Ref. f44 
argues against the presence of effects from soft radiation, their factorization formula does 
not have the analogue of the iSF. In our formalism the combined RG running of the two 
iBFs and the iSF cancels the running of the hard function as required by RG invariance. 
The presence of the iSF, which itself has a non-zero anomalous dimension, plays a crucial 



role in achieving this RG invariance as was shown in Ref. |17]]. Since Ref. [44| does not have 



the iSF they do not naturally achieve the required RG invariance. Instead, RG invariance is 
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implemented by introducing a 'hidden' Q 2 dependence in their two nucleon beam functions, 
which are individually ill-defined. This hidden Q 2 dependence is argued to arise from a 
collinear anomaly due to the absence of soft modes. It is further stated that soft modes 
have a vanishing contribution if the collinear anomaly is properly regularized. Since the 
presence of the collinear-anomaly already assumes the absence of soft modes, we do not find 
this argument compelling. The hidden Q 2 dependence reappears in their matching of the 
nucleon beam functions onto the PDFs at the p^-scale, despite the fact that p\ <C Q 2 . All 
of these problems are avoided if one starts with the correct degrees of freedom and includes 
the effects of both collinear and soft radiation. SCETn is known to be the appropriate EFT 
for this purpose. 



III. THE NON-PERTURBATIVE p T REGION 



In the previous section we reviewed the factorization formula for the region Aqcd 
Pt <C M. We now consider the pt distribution in the region where pt ~ Aqcd- The 
factorization theorem is given by 



d 2 a 7r 



dp 2 , dY N 2 



where !C q is defined as 



/ dx 1 dx 2 H q z (x 1 x 2 Q 2 , n Q ] (i T ) IC q (x 1 ,x 2 ,PT,Y, fi T ), 
Jo Jo 

(8) 



/f f d 2 b 
dt n J dt n J j^^Jo( b ±PT)B^(x 1} t^,b ± ,fiT)Bl(x2,t^,b ±} fi T ) 



x S~ q \ Xl Q - e Y Jp 2 T + M 2 - -±-,x 2 Q - e~ Y Jp\ + M 2 - -^,b ± ,/i T ). 

(9) 

In this case, the iBFs (B q n ) and the iSF (S^ 1 ) are evaluated at the scale \i T ~ p T ~ Aq C d 
with the hard function H q z evolved via its RG equations down to this same scale. Since 
fjbr ~ Aqcd, the iBFs and the iSF are non-perturbative. This expression for fC q was already 
derived in Ref. [H], |T5|] . In that work we focused on the region p T S> Aq C d so that the iBFs 
and iSF were perturbative and the iBFs were further matched onto PDFs. In this case, 
since fiT ~ Pt ~ Aqcd, the iBFs and iSF are non-perturbative. A perturbative matching 
onto PDFs is no longer valid and the final form of the factorization theorem is given by 
Eqs. (H) and @. For phenomenological purposes, the non-perturbative function K? must 
be modeled. When px ^> Aqcd, the scale \i? ~ pt is perturbative. K, q then becomes 
a perturbative object and the iBFs can be matched onto PDFs as in Eq. (0), leading to 
Eq. (D. 

We model the function K, q by imposing two requirements. First, the model for K, q must 
preserve the correct RG evolution properties so that it cancels the running of the hard 
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function H\ , as required by the scale invariance of the cross section. Second, as one increases 
p T from the non-perturbative region to higher perturbative values, Eq. (|8|) must reduce to 
Eq. (^). In order to smoothly transition between the non-perturbative and perturbative 
values oipx, we write the the iBFs in Eq. (|S]) as in Eq. (|^), even in the non-perturbative region 
where ~ Pt ~ ^qcd- In this region, Eq. ([?]) is no longer a perturbative matching equation 
but instead defines a new non-perturbative function X n . qr . As one increases /it ~ Pt to 
perturbative values, the function X n . qr corresponds to the perturbatively calculable coefficient 
in the matching of the iBF onto the PDF. Similar statements apply to the rz-sector iBF. 
With these conventions, one can write the function )C q as 

f 1 dx' f 1 dx' 

K. q (x l ,x 2 ,p T ,Y,n T ) = 22 — T~ l —rG qrs (x 1 ,x 2 ,x' 1 ,x' 2 ,p T ,Y,fi T )f r (x' 1 ,^ T )f s (x' 2 ,n T ), 

rs J xx x l Jx2 X 2 

(10) 

for all values of pt- For perturbative p^-values, the quantity Q qrs is perturbative and 
identical to that given in Eq. @. The expression in Eq. (|8|) then properly reduces to the 
factorization theorem of Eq. (|), valid in the region Aq C d < Pr <S M z . For fi T ~ p T ~ 
Aqcd, Q qrs is non-perturbative, as are all the quantities on the RHS of Eq. (||). In the non- 
perturbative-pr region, Q qrs can be interpreted as the non-perturbative TMF controlling the 
dynamics of transverse momentum dynamics of order Aqcd- The modeling of the function 
K. q in this region is reduced to the modeling of the TMF Q qrs . 

We view Eq. ( |I0D with K, q rewritten in terms of the standard initial state PDFs and a 
new TMF function Q qrs as more convenient than the form in Eq. (|9|). Both ways of writing 
IC q are equally valid. In Eq. (|]), the iBFs might be associated with TMDPDFs in the 
language used for the study of SIDIS processes. These iBFs are invariant under covariant 
gauge transformations but are not in general invariant under singular gauge transformations. 
However, the full product of the two iBFs and iSF that define JC q is completely gauge 
invariant. For a more detailed discussion of this point we refer the reader to Ref. [0. The 
form of K, q in Eq. (|l^) makes gauge invariance manifest. Since both K q and the PDFs are 
gauge invariant, the TMF Q qrs is also gauge independent. In Eq. (0), the gauge dependence 
of the iBF under singular gauge transformations is isolated into the function X n . qr . This 
situation also applies to the fi-sector iBF. The gauge dependence of the iBFs then cancels 
in the product that defines Q qrs . In this way, the non-perturbative dynamics in the region 
Pt ~ Aqcd is described in terms of gauge invariant initial state PDFs and the TMF. 

IV. TMF MODELS AND NUMERICAL RESULTS 

In this section we develop phenomenological models for the TMF function Q qrs in the 
non-perturbative region. We require that the model for Q qrs reduces to the perturbatively 
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calculable result as one increases px- We write Q qrs in the form 

g qrs ( y x 1 ,x 2 ,x' 1 ,x' 2 ,p T ,Y,fi T ) = dp' T g^ t (x u x 2 , x[,x' 2 ,p T ^l + (p' t /pt) 2 , Y, /x t ) 

x G mod (p' T , a, b, A), 

(11) 

which is a convolution of the partonic result for the TMF function Gp™t with a model 
function G mo d [£5-fT7fl. This form is reminiscent of that used in the CSS approach, where 
the integrand of the Fourier transform is decomposed according to 



W(b) = W(K)W NF (b), b 



V 1 + Wr, 



(12) 



W(b) is the perturbative resummed contribution and W NP denotes the non-perturbative 
contribution. b max is a free parameter typically taken to be of order 1 GeV -1 . 
We parametrize our non-perturbative contribution as 



G mod (p' T ,a,b,A) 
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A2 I a7 



a-1 



exp 



{Pt ~ 
2A 2 



(13) 



and fix N by the normalization condition 



dp' T G mod (p' T , a, b, A) 



1. 



(14) 



In principle, the model function G mo d can have flavor indices r, s. For the sake of simplicity 
we will work with a flavor-independent model function G mo d- Different choices of the pa- 
rameters a, k, A correspond to different model choices for the non-perturbative TMF Q qrs . 
The model function parameters are chosen such that G mo d will peak at p' T ~ Aq C d with 
an exponential fall off for larger values of p' T . As a result, Q qrs in Eq. (|ll]) receives sizeable 
contributions only from the region p' T ~ Aqcd- Thus, in the region px ^> Aqcd one can 
Taylor expand QZ^t. around the limit Pt ^ Pt ~ ^qcd- When combined with Eq. ( fU]) this 



gives 



G qrs (xi, x 2 , x[, x' 2 ,p T , Y, ht) 



^part.^l, X 2 , X[, X 2 ,p T , Y, fl T ) + 0( 



Pt ' 



(15) 



In the region of perturbative pt, the function Q qrs properly reduces to its perturbative 
limit with all model dependence suppressed by powers of Aq C d/pt- In this way, the model 
dependence is restricted to the non-perturbative region, as expected. The perturbative 
region of the pt spectrum remains calculable in a model-independent way to leading order in 
Aqcd/pt- One could consider more sophisticated model functions that contain ^-dependence 



S 




FIG. 1: The result for the py-spectrum of the Z-boson for the best fit parameter choices a 



2.25, b 

PTmin : 



0.1GeV,A = 0.45 GeV. We have also set y? Q 



-M| and fj, T 



PT + PTmin Where 



1 GeV. The data points were collected by the CDF and DO collaborations [22, 23 1. 



and that incorporate additional effects, but we restrict ourselves in this initial analysis to 
the form of Eq. (|l"3| ) . 

The implementation of the model also requires care regarding the choice of the scale 
\i T . In the perturbative pt region, the scale \it ~ Pt is the appropriate choice. However, 
one cannot use \xt ~ Pt when pt is of order Aqcd or smaller. The RG equations for the 
evolution of the hard function H q z (x\X2Q 2 , Hq] Ht) become non-perturbative in this region, 
and G^t. 111 Eq. ( |PT|) becomes incalculable. A sensible choice for \xt that can be applied in 
both the perturbative and non-perturbative pt regions is 

t4 = £ 2 PT + PTmin, ( 16 ) 

where PTmin ^ 1 GeV is a low, but still perturbative, scale and can be viewed as another 
parameter of the model. It is analogous to the parameter b max that appears in the CSS 
approach to transverse momentum resummation. £ is a scale variation parameter we take to 
be 0(1). The above choice of scale for /x^ has several useful properties. As pr — » 0, the scale 
Ht ~^ PTmin so that Q^art 111 Eq. flTTD is still evaluated at a perturbative scale. Similarly, the 
running of the hard function H q z {x\XiQ 2 , \iq\ \it) will freeze at the perturbative scale PTmin 
as pt —> 0. For larger values of px S> PTmin ^ 1 GeV, Ht ^ £,Pt so that the appropriate 
choice of \it ~ Pt in the perturbative region is recovered. 
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FIG. 2: The result of varying the model parameters a, b, and A within their 68% confidence level 
allowed region. We have chosen [Aq = — M|, \a\ = p\ + Pr-min W1 th prmin = 1 GeV. We see that 
the variation of the model parameters only affects the very low px region and has a negligible effect 
in the region px ^> Aqcd- The data points are from the CDF and DO collaborations [22, B3J. 



We now present an example fit of the TMF function Q qrs to Tevatron data for the Z- 
boson p T spectrum. We choose [Aq = — Mf [|8| ^5[| , il t as in Eq. fllBD with £ = 1, and 
for simplicity set prmin = 1 GeV. We note that this ensures that the scale ht at which the 
PDFs are evaluated always remains at or above the initial scale Qo — 1 GeV used in the 
MSTW fit |50|| , a criterion pointed out in previous work in the CSS approach j3T|. We then 
perform a chi-squared fit of the parameters a, b, and A in Eq. (|T3|) against CDF data [52]; for 
simplicity we do not include the DO data in this example fit. The best fit values obtained 
are a = 2.25, b = 0.1 GeV, A = 0.45 GeV with a goodness-of-fit measure ^/d.o.f. ~ 0.7. 
The result for these best fit values are shown in Fig. [I] along with the CDF and DO data 
points. Fig. [l] shows that the TMF model is flexible enough to give a good description of 
data in the region p? < 1 GeV where non-perturbative transverse momentum dynamics 
becomes important. At the same time, a good description of the data is also achieved 
for larger perturbative values of pt where the result is given in terms of a perturbatively 
calculable TMF function. The model dependence introduced by G mo d turns off in the region 
Pt 3> Aqcd, as expected. This is further illustrated in Fig. |] where we show the results for 
the 68% confidence level region in the parameters a, b, and A. We see in Fig. [2] that while the 
different parameter choices affect the pT-distribution in the non-perturbative region, there is 
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almost no effect in the region 3> Aqcd- This is a reflection of Eq. fll5|) which shows that 
for p T ^> Aqcd the model dependence is power suppressed and the TMF function reduces 
to the expected partonic result. 

Before concluding we comment briefly on the universality of G mo d- We have neglected 
the possible flavor dependence of this function in our fit, indicating that we expect the non- 
perturbative dynamics of the valence up and down quarks that dominate Z-boson production 
at the Tevatron to be the same. However, it remains to be seen whether the valence-sea 
scattering which occurs in pp collisions can be described by the same G mo d- For this reason 
we refrain from making predictions for LHC production until this universality is tested by a 
detailed fit to the available data. We note that W NP in the CSS approach has been found 



to satisfy the universality assumption [31 



V. CONCLUSIONS 



In this manuscript we have performed an initial analysis of the Z-boson transverse mo- 
mentum distribution in the region p? ~ Aqcd using a factorization and resummation the- 
orem derived in SCET. Combined with our previous work |~7|, [TJ§ which focused on the 



region Aq C d <C Pt ^ M, a description of the entire p^-spectrum is now achieved in the 
framework of SCET. This formalism is free of the Landau poles that arise in the traditional 
approach to low-pr resummation in impact-parameter space, and are therefore independent 
of ambiguities and numerical difficulties which arise when transforming back to momentum 
space. In the region where px ~ Aqcd, the transverse momentum spectrum is affected by 
new non-perturbative effects that cannot be described by the standard PDFs alone. A new 
transverse momentum function (TMF), fully gauge invariant and defined in SCET, arises in 
addition to the standard PDFs. The TMF captures the non-perturbative dynamics associ- 
ated with the initial state transverse momentum distributions and with final-state emissions 
having transverse momenta of order Aqcd- We have devised phenomeno logical models for 
the TMFs in the region ~ Aqcd- These models are such that the TMF reduces to the 
expected perturbative result when pt ^> Aqcd- This allows for a smooth transition between 
the non-perturbative and perturbative values of pt- The TMF models also have the correct 
renormalization group evolution properties built in. We have given example fits of the TMF 
model to Tevatron data. The results of the fit for the TMF function give a good description 
of the CDF and DO data over the entire px spectrum. 

The work presented here is simply the first step in understanding the non-perturbative 
transverse momentum region within SCET. A more global analysis of the available data is 
left to future work, as is the modeling of the TMF for pp initial states. In principle, the 
TMF is different for pp and pp initial states. The universality of this function remains to be 
studied. These questions must be addressed to present predictions for the pt distribution 
at the LHC. We look forward to these future investigations. 
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